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We performed an exact diagonalization study of the spin 
liquid phase of the Multiple-Spin Exchange model on the tri- 
angular lattice. It is characterized by no Neel Long Range 
Order (NLRO), short- ranged magnetic correlations and a spin 
gap. We found no LRO in any local order parameter we inves- 
tigated (chiral, dimer,...). The probable asymptotic ground 
state degeneracy is discussed. We argue that it could be of 
topological origin and that the system is probably not a chiral 
spin liquid. A possible relation to the Afneck-Kennedy-Lieb- 
Tasaki (AKLT) phase is discussed. 

PACS numbers: 75.10.Jm; 75.50.Ee; 75.40.-s; 75.70.Ak 



I. INTRODUCTION 

A general model. The Multiple Sain Exchange 
(MSE) model was introduced by ThoulessH It is an ef- 
fective description of the magnetic properties of quasi- 
localized fermions. This lattice spin Hamiltonian was 
introduced to describe the magnetic consequences of tun- 
neling events where several indistinguishable particles ex- 
change their positions. For spin-^ particles, when two- 
body processes are the only important events, it leads to 
the familiar Heisenberg Hamiltonian: 

h = E -M'..- = E J ^ ( 25? > • $i + 1) (!) 

where P^j is the spin permutation operator, and 2Jjj the 
(positive) frequency of the tunnel process exchanging the 
particle on site i and the particle on site j. But in low 
density crystals, cyclic exchange processes P n of n > 3 
particles can be important. It is the case in low den- 
sity solid 3 He 64ms (□ and references therein) and in the 
Wigner crystafl. More generally, multiple-spin interac- 
tions with n > 2 are expected in the spin systems where 
interactions are strong and cannot be reduced to Heisen- 
berg couplings. In numerous quantum crystals, such cou- 
plings are expected to be generated by integrating out 
the non-magnetic degrees of freedom. For instance, the 
role played by MSE in the context of metal-insulator 
transitions has been discussed recently by Chakravarty 
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et alB. Therefore it is of wide interest to understand the 
nature of the magnetic fluctuations introduced by MSE 
couplings. 

The MSE spin liquid phase. — In a previous work we 
studied the phase diagram of the MSE model on the tri- 
angular lattice and showed the existence of a spin liquid 
phaseQ. In this work we present a detailed characteri- 
zation of the low energy physics of this new spin liquid. 
The hamiltonian involves the five simplest ring exchange 
patterns on the triangular lattice: 

h = j 2 J2p2-J 3 E ( p 3 + ^ 1 ) 

~ A 

+j 4 E (Pi+p* 1 )-^ E ( p 5 + p 5 _1 ) (2) 

+j 6 E (Pz + Ps 1 ) 

O 

Due to the Pauli principle, exchange of an odd number 
of fermions is ferromagnetic whereas exchange of an even 
number if antiferromagnetic. This is the reason for the 
alternating signs in Eq. (0) (all the J n being > 0). Be- 
cause of the equality P ia , 3 + ^3,2,1 = Pi,2 + P-2,3 + P3,i-^, 
valid for spin-^, triple exchange around triangles can be 
taken in account just by modifying the bare J2 frequency 
J 2 — > J2 — 2J3. In the following we thus assume without 
any loss of generality that J3 = and that J2 can be > 
or < 0. 

At T = there is a first order phase transition be- 
tween a paramagnetic phase (the ground state is a singlet 
S = 0) and a ferromagnetic one (the ground state is fully 
polarized: S — N/2)a. In this work we are interested 
in the ground state and first excited states structure in 
the paramagnetic phase. The behavior of this spin liq- 
uid when an external magnetic field is applied will be 
discussed elsewhere. 

What kind of Spin Liquid is the MSE model ? — 
Many different quantum ground states, which differ by 
their spatial broken symmetries and excitation spectra, 
have short-ranged spin-spin correlations and an S = 
ground state. Among the two-dimensional Hamiltonians 
with a single spin in the unit cell (no dimerization nor 
inequivalent bonds) there are three important examples 
where the ground state wave-function is understood: 

• Spin— I Klemi models. This family- of Hamiltoni- 
ans generalizes the Majumdar-GosliD model. Since 
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any first-neighbor valence-bond state is a ground 
state, the degeneracy is extensive in two dimen- 
sions. 

• Valence Bond Crystal: It has some dimer-dimer 
long-ranged order or some more complicated pla- 
quette order. It has a spin gap and spin-spin cor- 
relations are short-ranged. The translation sym- 
metry is spontaneously broken and therefore the 
ground state is degenerate. An example may be [the 
frustrated antiferromagnet on the square latticeS'tl. 
These spatial symmetry breakings-, a*e found in 
large-TV limits of SU(N) modelsMlJ when the 
spin 5 does not match the coordination number 
z (25 ^ mod z). 

• Valence Bond Solid (Affleck et alM). It exists 
when the spin 5 on one site matches the z coor- 
dination number of the lattice (25 = z). It breaks 
no symmetry and has a spin gap but, a priori, it 
cannat_.be constructed with spin-i. The Haldane 
phasdij in one dimension (integer spin) belongs to 
the same universality class. The ground state is 
non-degenerate (infinite volume limit with periodic 
boundary conditions), spin-spin and dimer-dimer 
correlations decay exponentially with distance. 

There are very few models where resonances between 
short-ranged valence-bond states select a single (and 
fully symmetric under symmetry operations) combina- 
tion of them_A perturbed Klein Hamiltonian could be 
a realizationtJ. The J1-J2-J3 model on the square lat- 
tice might also be.-a. candidate (exact diagonalizations on 
a N = 16 sampleE!.). So far, the short-ranged Resonat-. 
ing Valence-Bond (RVB) picture proposed by AndersonEZl 
has not yet found an explicit realization: there is no def- 
inite spin-i Hamiltonian in two dimensions with no bro- 
ken translation symmetry and a non degenerate ground 
state. From this point of view the MSE is of great rel- 
evance since our numerical data point to short-ranged 
correlations with no kind of LRO. 



II. RANGE OF THE SPIN LIQUID PHASE IN 
THE MSE MODEL ON THE TRIANGULAR 
LATTICE 

Even in the classical limit, the ground state of the MSE 
Hamiltonian is exactly known only in three limits: 

• pure J2 > case: 3-sublattice Neel State 

• pure J4 case: 4-sublattice Neel State i±e±rahedral 
state found by Momoi, Kubo and NikO'Ea). 

• pure J n with odd n: Ferromagnetic. 

An approximate phase diagram for the classical system 
at T — has baea obtained in the variational subspace 
of planar helicesE3, and in the 4-sublattice subspace for 



the J2 — J4 model (Kubo et a/Jia). From these results we 
can sketch the qualitative shape of the classical diagram 
(dotted line of Fig@). 

In the 5 = „ quantum system there is hardly any exact 
result, but some regions of the diagram are understood: 

• Pure J2 > case: the A.F.Jieisenberg Hamiltonian 
has 3-sublattice Neel RLOB. 

• Large region about the pure J4 case: no Long 
Range Order (LRO), finite -spin gap and short- 
ranged spin-spin correlationsB. 

• Large ferromagnetic phase including the pure J2 < 
and pure J5 Hamiltonians. 

A possibility for the quantum phase diagram is pre- 
sented Fig. |l|. This simple guess relies on the following 
data: 

• Preliminary exact diagonalizations results indicate 
that the extension of 3-sublattice Neel phase in the 
J2 — J4 model is strongly reduced by 4-spin ex- 
change. This is already-jthe case at the mean field 
Schwinger Boson leveled. The possibility that this 
Neel phase is destroyed by an infinitely small J 4 is 
not excluded (additional work is in progress). 

• Our exact diagonalization data show no LRO in the 
large- J4 region. 

• The energy of the ferromagnetic state is the same 
in the classical and quantum cases but A.F. con- 
figurations gain energy from quantum fluctuations. 
Therefore, the ferromagnetic phase is slightly re- 
duced in the 5 = | case. 

On the J2(> 0) - J5 line, we suggest the existence of 
a spin liquid window between the ferromagnetic and the 
3-sublattice Neel phase. However, this has not yet been 
numerically investigated. Another possibility is a first or- 
der phase transition between NLRO and ferromagnetism. 
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FIG. 1. Qualitative phase diagram of the MSE model on 
the triangular lattice. The right corner is the pure J4 model, 
the left one is the pure J5 one. Down and up corners are the 
ferromagnetic and antiferromagnetic Heisenberg Hamiltoni- 
ans. Letters and dashed lines: Schematic phase diagram 
for the classical model, a: Ferromagnetic, b: 3-sublattice Neel 
state, c: 4-sublattice Neel state, d: long wave length spiral 
state, or ferrimagnetic state. Full lines: likely scenario for 
the S = I quantum phase. We put a small disordered window 
on the upper J4 = line but a first order transition between 
ferromagnetism and 3-sublattice Neel order is also possible, 
some more numerical work is needed to fix this point. 



III. SHORT-RANGED MAGNETIC 
CORRELATIONS 

A. no Neel Long Range Order (NLRO) 

The first question to address in the non-ferromagnetic 
region is: Is the system Neel Long Range Ordered at 
T = ? 

Periodic boundary conditions. - As a Neel state 
breaks the SU(2) and some spatial symmetries of H, it 
cannot be an exact eigenstate on a finite size system. It 
has to be a linear combination of eigenstates belonging to 
different Irreducible Representations (IR) <j*f-ifte spatial 
symmetry group and to different S sectorsEilES. As the 
dynamics of a Neel order parameter is the one of a free 
rotator, the corresponding low energy levels scale as: 



E 



S(S + 1) 
Nxo 



(3) 



since the inertia of that rotator is proportional to the 
number of sites N (\o is the susceptibility per site at 
zero field: xa — jj 9< aB > 1 r-q)' ^ nxe d S, these states 



N dB \B = 

collapse to the ground state in thermodynamic limit more 
rapidly that the softest magnor£J. These states form a 
tower of states. This low energy structure is absent in 
the spectra as soon as J 4 is not negligible. Fig. HJ1 shows 
how the 3-sublattice Neel tower of state is destroyed by 
4-spin exchange. 



The softest magnon has a total energy E„ 



= c|k„ 



where c is the slowest spin-wave velocity and k m i n ~ -h= is 
the smallest wave vector allowed by the periodic boundary 
conditions on the finite size sample. This gives E m in — -^=f- 
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FIG. 2. First-neighbor Heisenberg spectrum on the trian- 
gular lattice. The SU(2) and spatial symmetry breakings due 
to the 3-sublattice Long Range Neel Order appear as a set of 
states with an energy scaling as E(S) ~ S(S + 1)/N (dashed 
line). The symbols represent the quantum number of the dif- 
ferent eigenstates. k is the wave vector, k = B is the corner 
of the Brillouin zone (the impulsion is not indicated if k 7^ 
and 7^ B). Rg is the phase factor obtained in a ^-rotation 
about the origin and a stands for a reflection about an axis. 
The vertical lines indicate the energy range where the density 
of states is high and all the eigenvalues have not been com- 
puted. This has no consequence on the low energy part of 
each irreducible representation of the symmetry group, where 
the eigenstates are known exactly. 
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FIG. 3. Heisenberg spectrum perturbed with 4-body ex- 
change J4/J2 = 0.1. The Neel structure is destroyed: the 
states formerly embedding the symmetry breakings in the 
J4 = case have been "pushed up" in energy to the con- 
tinuum of excitations. 



Twisted boundary conditions. — For the tower of states 
to appear on a finite size sample, the sublattice structure 
must be compatible with the boundary conditions. An 
order with a long wavelength helix, or an incommensu- 
rate phase is therefore difficult to detect on small samples 
with periodic boundary conditions. Fortunately, twisted 
boundary conditions^ allow to overcome this difficultyElJ. 
For N — 19 and J2 = —2, J4 = 1, we scanned the Bril- 
louin zone to determine the twist qo which minimizes the 
ground state energy. The twist vector qo lies inside the 
Brillouin zone and has no particular symmetry. If the 
system had some NLRO (commensurate or not), this qo 
would indicate the ordering wave vector and the spec- 
trum would show a tower of states E ~ n%1 ■ Fig- % 
shows that it is not the case. 

N-19 J|"=-2 J 4 =l 
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FIG. 4. Spectrum with the twisted boundary conditions 
which minimize the ground state energy. Eigenstates repre- 
sented by black triangles are k = states. The spectrum 
shows no long-ranged (spiral) order: if there was a tower of 
states the lowest k = states (joined by the dashed line) 
would be lower than the other excitations and their energies 
would be proportional to S% . Notice that with these twisted 
boundary conditions the spatial symmetries (except transla- 
tions) as well as the SU (2) symmetry are lost. But S z is still a 
conserved quantity. The twist vector is qo = 0.27Ai— O.2OA2, 
where Ai and A2 are middle of the boundary of the Brillouin 
zone. By symmetry, there are twelve twists equivalent to qo. 

B. Spin-spin correlations: < St ■ Sj > 

For spin-7p Si • Sj varies between — 3/4 (singlet) and 
1/4 (triplet) and the statistical average is zero. Com- 
pared to these extremal values, the spin-spin correlations 
measured in the ground state of J2 = — 2 J4 = 1 for 
N = 16,24 and 28 are small (Table @). Fig. | displays 
the absolute value of the spin-spin correlation as a func- 
tion of distance. On these three samples the available 
distances are rather small — j\ < 3) and the data show 
important size and orientation dependence (for N = 24, 
the system has no K21 symmetry and there are 3 non- 
equivalent directions) . 
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TABLE I. Correlations < S t ■ Sj > in the MSE ground 
state of N = 16, 24 and 28 samples (J 2 = -2,J 4 = 1). These 
data are plotted in Fig. |E| The TV = 24 sample is a 6 x 4 
lattice which has not the IZ^ symmetry. The three directions 

are thus non equivalent, this explains why < S, ■ Sj > has 
three possible values at distance 1. The underlined values 
are over-correlated due to periodic boundary conditions: they 
correspond to antipodal sites on the torus. 



2 With periodic boundary condition, the spin Si at site i 
interacts with the spin Sj at site j. Using the wave vector q 
to twist the system makes the rotated spin TZ^.i interact 

with IZq.j j^Sj j . IZg is the rotation of angle 6 about the quan- 
tification axis z. Such an Hamiltonian is no longer SU{2) in- 
variant when q / 0: §t ot — §ij is no longer a conserved 
quantity but Sl ot commutes with the twisted Hamiltonian. 
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more regular and the rapid decay suggests a rather short 
correlation length. 

The point is now to understand the local structure of 
the ground state wave function. There are two naive 
ways to build an S = state out of a large number 
of spins- i: 1) combine ferromagnetically the spins in a 
small number of sublattices (NLRO) and let these large 
spins screen themselves. 2) combine the spins two by two 
in singlets. The absence of tower of states and significant 
ferromagnetic correlations have excluded the first possi- 
bility (no NLRO). It is easy to check that the weakness of 
the A.F. correlations makes the second scenario unlikely: 
the screening of a single spin at the origin involves an im- 
portant number of neighbors, up to distance d ~ 2 • ■ ■ 3 
(This is easily checked from the data of Tab. ||). This will 
be confirmed in the next paragraph where we show that 
dimer-dimer correlations are weak. 



C. Dimer-Dimer correlations 



FIG. 5. Absolute value of the spin-spin correlation (data 
of Table ^). The vertical range goes from to |. — | is the 
first-neighbor correlation in the Majumdar-Gosh chain. This 
scales has been chosen to show how small the correlations are 
in the MSE system. The filled symbols correspond to the 
correlations for which finite size effects are expected to be the 
smallest (see text). 

At this point, it is difficult to draw conclusions on the 
spin-spin correlation decay from Fig.El However, a finite 
size system with periodic boundary conditions is a torus 
and we expect some geometrical effects for pairs 
of sites which have particular symmetries. The strongest 
effect is the enhancement of < Si • Sj > on samples where 
j i and i j are equivalent vectors: antipodal sites 
are over correlated (underlined values in Table |). If the 
sample has no IZ^ symmetry the finite size effect on 



< Si ■ Si > at distance d 



3\ 



should be weaker 
in the direction of the vector v = i — j which is the not 
frustrated by the periodicity of the torus|. If we eliminate 
antipodal sites and frustrated directions in N = 24 only 
the filled symbols of Fig. [sj remain. The behavior looks 



A simple way to evaluate how the the v step is "frus- 
trated" by the periodic boundary conditions of the torus is 
to consider the closed loop starting from the origin: 0,0 + 
v, O + 2v, • • • , O + (n — l)v, O. We associate frustration with 
the winding number of this closed curve. The minimum num- 
ber is one: by successive translation of v we are back to the 
origin after winding one time around the torus. We make the 
assumption that unphysical effects due to finite size of the 
torus are smaller in the direction where the winding number 
is the smallest. Given the distance d, the relevant < Si ■ Sj > 
is thus obtained in the direction which gives the topologically 
simplest loop. 



We define the dimer operator on a pair of sites 
by d,., = 



This projector is 1 on a singlet and 



on a triplet. The dimer correlation between a reference 
bond (1,2) and is Du =< #|di )2 di,j|* > - < 

\I>|di, 2 |\I> >< 3>|c2j,j|* >. Q As for the normalization, 
we look for the maximum value of D^j. It is achieved 
when the two bonds are completely correlated and gives 
D itj =< *|di, 2 |* > - < *Kj|* >< I'Mi^l* >• So, 
we measure dimer correlations by: 

„, . = I hA (a) 

Pl '° < *|d 1)2 |* > (1- < tfldyl* >) 1 1 

_ < jf?\di, 2 di,j\% > - < vEKjjjE >< gjdi^aTg > 
(l-<*|d M |*>)<*|di, a |*> 

Pij is represented in Fig. ^| and |?]. For this quantity, 
zero means that the presence of a singlet on (1,2) and 
the presence of a singlet on are uncorrelated. A 

value of one means that if a singlet exists on (1, 2) there 
is always one on The minimal possible value is 

<di - J> ic „min _ _ .469 at distance 1 
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in the N = 28 ground state. We observe negative values 
on the 4 bonds which are at distance 1 from the refer- 
ence bond. Compared to p™ m , this values are important 
(30% ) and increase with the system size. This testifies an 
important short-ranged repulsion of parallel dimers and 
excludes the possibility of spin-Peierls like structures in 
which the lattice would be regularly tiled with parallel 
dimers. Fig. ^ shows that the local dimer distribution fa- 
vors angle of ±^ between dimers bonds. This tendency 



4 If the system is invariant under rotation of 27r/3 < 
$ | di,j \ty > only depends on the bond length \j — i\. 
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to a ±? ordering of dimers already appears in the 6- 
site system with open boundary conditions. This sys- 
tem (shape of a triangle) is the smallest system with a 
low ground state energy and a significant spin gap. In 
larger systems, this local ±-| structure is all the more 
pronounced as the ground state energy is low. 



FIG. 7. Values of pij on the N = 28 lattice: the reference 
bond is (1, 28) and the 4 most strongly correlated bonds are in 
full line, they are at distance d = \fl /2 from (1, 28) and form 
a triangular pattern. We notice that the 4 "first-neighbors" 
of the bond (1, 28) are strongly anti-correlated (see values on 
(6,10) (10,2) (23,27) and (18,23)). 
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distance 

FIG. 6. Dimer-dimer correlations in the ground state \^ > 
of the MSE exchange Hamiltonian in the Spin-Liquid phase 
(J2 = —2 and J4 = 1). They are plotted as a function of 
distance between bonds for N = 16 and 28. 



From Fig. ^ one could doubt that the dimer-dimer goes 
to zero with distance. However, we think it is the case 
for the following reasons: 1) dimer-dimer correlations are 
very weak. Even at short distances, the wave function 
has to be seen as a gaz or a liquid and not as solid of 
dimers. 2) The correlation at distances d > 2 decreases 
in a significant way from N — 24 to N = 28. 3) It 
is not possible to tile the triangular lattice with dimers 
at ±^ from their neighbors without defaults (one site 
out of seven would not touch any dimer, see Fig. ||) and 
this local property cannot be propagated to the entire 
lattice. We stress that this ±^ dimer correlation cannot 
picture the entire wave function. It would not account 
for the low first-neighbor "dimer density" : < Si- Sj > at 
distance \i — j\ = 1 is only very weakly antifcrromagnetic 
(< S t ■ Sj >~ -0.07) and for N = 28, < d h2 > is even 
lower than the second-neighbor one < d\ t \\ >, see Table 
[j].). The information provided by the dimer structure 
of Fig. is that a low energy state cannot be achieved 
without a local order with a complex geometry. 
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FIG. 8. Dimer covering of a 6/7-depleted triangular lattice. 
It is not possible to tile the triangular lattice with dimers 
forming angles ±7r/3 without defects. On site out of 7 does 
not have any dimer (circles). It is worth noticing that a multi- 
ple-spin exchange Hamiltonian for which this state is an exact 
ground state can be constructed by the procedure introduced 
by KleinQ. The Klein Hamiltonian is a sum of projectors: 
each projector acts on a rhombus and projects the four sites 
in their S = 2 subspace. All the sites are equivalent on the 
depleted lattice and the Hamiltonian can be expressed as a 
Multiple-Spin Exchange Hamiltonian containing Pij at dis- 
tance 1 (J2 = 5) and %/3 (JJ 11 = 1) and 4-body terms on 
rhombus ( J4 = 1 for Pi, 2, 3, 4 4- h.c. and J 4 = 2 for a Pi, 3^2,4 
term on the same rhombus). Actually, this state is also a 
ground state for the first-neighbor antiferromagnetic Heisen- 
berg Hamiltonian on this lattice. 
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D. Chiral-chiral correlations: < Im(P; j i k)Im(Pii ji t y) > 

Kubo et al. have demonstrated that the classical 
ground state of the J4 — 1 Hamiltonian was a 4-sublatticc 
Neel state, where the magnetization Abactors of the 4 sub- 
lattices form a regular tetrahedronl 18 H 19 l. They showed 
with Monte-Carlo simulations that the classical-system 
has a chiral phase transition at finite temper atureij. This 
transition is associated with the ferromagnetic Ising-likc 
order parameter defined on three sites by the operator: 
k = 2(Si x S3) ■ §2- 111 the classical ground state all tri- 
angles have a chirality k = or k = -y. Kubo's 
simulations indicate that the < t(i,2, 3)^(1' ,2',3') > corre- 
lation function remains long-ranged up to a finite critical 
temperature. To check this hypothesis in the quantum 
case, we have computed this correlation in the ground 
state of the quantum system for the pure J4 hamiltonian 
(like Kubo et al.) and found a significant decay with 
distance (Fig. HID). Moreover, this correlation shows 
negative values, which is unexpected for a ferromagnetic 
Ising model. The existence of such Ising LRO in the 
quantum system, even at T=0, seems unlikely. 



expectation value on 3 free spins: Tr[K 2 ]/Tr[l] = 0.375. 
On the other hand, it is both much smaller than the 
maximum quantum value of n 2 (= §) and much smaller 

aw lying singlets on the 
3. The last reason why 



than the value obtained in the 
Kagome lattice (< n 2 >~ 0.7)1 

the value of < k 2 > should not be interpreted as sufficient 
evidence of T-symmetry breaking is that it contains no 
information but the first-neighbor correlation: 



< k >= 



< Si ■ Si > 



(5) 



We also performed computations at the J2 = — 2 and 
J4 = 1 point. Unsurprisingly, the ferromagnetic coupling 
reduces the chirality. In the ground state of N — 28 at 
J 4 = 1 and J 2 = -2 we found < k 2 >= 0.4791. This 
is readily understood from Eq. B. Since < Si ■ Sj > 
is enhanced by a negative J2, local fluctuations of chi- 
rality are reduced. It is more interesting to look at 
< Im[P]Im[P'] > as a function of distance. The data 
are shown in Tab. || and support an absence of long- 
ranged chiral order at the J 2 = — 2 and J4 = 1 point as 
at the pure J4 point. 
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FIG. 9. Chiral-chiral correlation function for the J4 = 1 
Hamiltonian. The distance is measured in lattice spacing 
units between the center of gravity of the two triangles. The 
chirality is simply expressed with spin exchange operator: 
lm(P iJt k) = 5j (P< j,* - Pk,j,i) = 20* x S k ) ■ Sj = k. The 
value in the classical 

tetrahedral state is < Im[P]Im[P'] >= i ~ 0.037 for any 
couple of clockwise triangles P and P' (dotted line). If the 
two triangles have some sites in common, < Im[P]Im[P'] > 
may have a (small) imaginary part. In this case, the real part 
is plotted. 



|(123) 


- (1'2'3')| 


J 2 = -2 J 4 = 1 


J 2 = J 4 = 1 




1 


-0.0285 


-0.0306 


3 


= 0.527 


+0.0280 


+0.0356 




2 


-0.0078 


-0.0132 


V7 


= 2.646 


+0.0018 


+0.0009 



TABLE II. Chiral-chiral correlation in the ground state 
of an N = 28 system. < Im(Pi i 2,3)Im(P 1 /. 2 '. 3 ') > at J 4 = 1 
and for two values of J2. (1,2,3) and (l'2'3') are clockwise 
triangles and | ( 123) — (l'2'3')| is the distance between the 
centers of the two triangles. The increase of the chirality 
at distance y/7 is very likely to be a finite size effect. For 
the three other distances we observe a significant decrease 
of chirality due to the first-neighbor ferromagnetic coupling 
(Ja = -2). 



It is true that the local fluctuations of k are large: for 
N = 28 (T = 0) < n 2 >= 0.5221 which is larger that the 
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IV. LOW ENERGY SPECTRAL PROPERTIES 

In the previous section we have reported some correla- 
tion data showing that the ground state wave function of 
the MSE Hamiltonian in the spin-liquid phase has short- 
ranged correlation functions. We now describe the low 
energy spectrum of the system. We emphasize on the 
spin gap and the broken symmetries. Again, we concen- 
trate on the J2 = — 2, J4 = 1 Hamiltonian. 

A. Spin gap 

1. Finite size effects on the ground state energy 

For small sizes the spin gap and ground state energy 
are sensitive to the sample shape and size. For this rea- 
son, an N = 00 extrapolation of the spin gap value 
is not straightforward and requires a precise analysis. 
The apparently irregular behavior of the spin gap and 
ground state energy for systems with less than N ~ 24 
sites (Fig. [lO]) indicates the short-ranged rigidity of the 
ground state wave function. This agrees with the pre- 
vious section: the characteristic length £ of the ground 
state correlations is a few lattice sites and for systems 
of this size, the spectrum is sensitive to boundary condi- 
tions. This effect can be used to probe the local struc- 
ture of the ground state wave function and to identify 
some short wave length resonances which lower its en- 
ergy. Two families of samples should be distinguished: 
systems with even N have a low ground state energy and 
a higher spin gap that the systems with N odd. Within 
each family, the behavior with the size is more or less 
monotonous: for even-N systems (resp. odd-N) the en- 
ergy per site increases (resp. decreases) with N and the 
spin gap decreases (resp. increases) with N. The impor- 
tance of a 6-site periodicity for achieving a low energy^ 
agrees with the ±7r/3 orientation found by looking at 
the short-ranged dimer-dimer correlations (Fig. Q). The 
enhancement of this ±-| structure in the lowest energy 
systems shows that it captures some important short dis- 
tance correlations. 

An important point is that both families, odd and 
even, merge for sizes N > No ~ 30, (either if we look 



at E/N or at the spin gap). We have interpreted this 
feature as a crossover behavior: above the characteristic 
size Nq ~ £ 2 the finite size correction decay exponen- 
tially fastQ. Fig. [[2] also shows a rapid decay of finite size 
effects when approaching N — 36. From this argument, 
the N = 36 system can almost be considered as having 
reached the thermodynamic limit. This is all the more 
probable as this sample (6 x 6) has the symmetries of the 
infinite lattice and does not frustrate any short-ranged 
order that we observed to be favorable. 



2. Spin gap 

The correlation between the ground state energy and 
the energy gap to the first S = 1 (and S = § , if N is odd) 
excited state is indeed precise (see Fig. [ll]) : the most sta- 
ble systems (i.e. with E/N minimum) are those with the 
largest spin gaps. The different spectra roughly fall onto 
one line when the spin gap is plotted versus E/N. This 
suggests that the position on this line is a single param- 
eter which measures the frustration of a system. On the 
right part of Fig. [ll], systems have important frustra- 
tion, because of geometric constraints the local structure 
of the wave function which would minimize the energy 
in the infinite systems cannot develop. Systems on the 
left part benefit of particular boundary conditions which 
allow some stabilizing resonances and therefore an E/N 
lower than in the thermodynamic limit. From the data of 
E/N versus A we can estimate the position of the infinite 
system. E/N converges to —3.95 ± 0.05 (safe estimate 
from Fig. [H]). With this energy interval, the value of the 
spin gap can then be extracted from Fig.[ll| We obtain 
A = 1.3 ±0.5. 



6 Notice that the situation would be different in the case of 
NLRO. The finite size corrections of the ground state energy 
per site are proportional to 7V -3//2 (linear spin-wave result). 
The energies of the different families of samples (those which 
frustrate the LRO and those which do not) would intercept 
linearly at N = 00 in a plot of E/N versus N' 3/2 . As for the 
spin gap, it would vanish as N~ : this is definitely inconsis- 
tent with present results (Fig. fid) 



5 Systems with the lowest ground state energy per site turn 
out to be those with wave vectors k = A, B and B/2 in their 
first Brillouin zone. An illustration of this sensivity to the 
shape is provided by the two N = 24 systems of Fig. [ll]. The 
more stable N — 24 system is the 6x4 torus (which possesses 
the four wave vectors mentioned and E/N — —4.03). It can 
be compared with the most compact 24-site torus (label 24' in 
Fig. [ll] and empty square at N = 24 in Fig. [lo|). The later has 
only one k = A vector and is more frustrated: E/N = —3.88. 
The same behavior was observed for 12 sites. 
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B. Ground state symmetries in the thermodynamic 

limit 

The spin gap described in the previous paragraph in- 
sures that in this region of the phase diagram the system 
does not spontaneously break the SU(2) symmetry at 
low temperatures. We now turn to the 5 = sector and 
address the question of the ground state degeneracy and 
the possibility for a spatial group symmetry breaking. 
This is important to characterize the class of spin liquid 
this model belongs to. 



1. Singlet states in the spin gap 

In nearly all even samples (TV =36, 30, 28, 24, 18, 16, 
12) the ground state belongs to the trivial representation 
of the point group and translation group. This represen- 
tation being one-dimensional, the ground state is non- 
degenerate for finite TV. Above the ground state we find 
always less than 10 eigenstates in the spin gap (i.e. below 
the first triplet energy). There is no extensive entropy at 
T = nor any singlet soft mode (notice the difference 
with the Kagome lattice Heisenberg antiferromagneto) . 
This implies that the system at J 2 = — 2 and J4 = 1 is 
not a quantum critical system. 

We have analyzed the (spatial) quantum numbers of 
the low energy singlet states to detect if some symmetry 
sectors were collapsing to the absolute ground state. So 
far, because of the poor regularity in the symmetries, this 
task has only brought partial results. Here are the most 
important data to be understood: 

• There is no low-energy singlet in the k = B 
sector (±B are the two corners of the Brillouin 
zone). This excludes any 3-step translation sym- 
metry breaking. It also confirms that the ±7t/3 
dimer order (Fig. 0), which looks 3-step periodic, 
is only local. 

• The frequent occurrence of k = A states in the 
spin gap suggests that such states collapse to the 
ground state and give a 2-step translation symme- 
try breaking (A is one middle of the boundary of 
the Brillouin zone). In particular the energy gap 
between the absolute ground state and the first 
k = A state drops by a factor of ten between 
TV = 24 and TV = 36, (2.5 at TV = 24 and 0.243 
for TV = 36 see Pig ph. 



• The TV = 36 singlet sector of the J2 = — 2 and J4 = 
1 Hamiltonian has a quasi four-fold degeneracy of 
the ground state (Fig. |l3|). The lowest energy is 
a single state and the next energy level, which is 
immediately above, is 3-fold degenerate. Because 
of the large size of the TV = 36 system (compared 
to the supposed correlation length), we believe that 
this feature is a particularly valuable information 
on the spectrum of the infinite system. 
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In the next paragraphs we analyze these results. First, 
we argue that these data are not in favor of "spin-Peierls" 
like symmetry breaking. Then we show that a four-fold 
degeneracy could be of topological origin. Thirdly, we 
explain why the system is probably not a Chiral Spin 
Liquid. Eventually, we discuss the link between the MSE 
spin liquid and some kind of AKLT (or Haldane) phase 
on a square superlattice with four spins in the unit cell. 
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glet excited state (crosses). 



3. Topological degeneracy ? 

But this degeneracy matches the prediction of a topo- 
logical degeneracy due to the periodic boundary condi- 
tions and to the non-trivial topology of the torus. 

Most of the important arguments for this topological 
degeneracy were proposed by Rokhsar apd KivelsonEJ, 
Read and Chakrabortyc3, and SutherlandcS but here we 
sketch a more formal demonstration. 

Dimer and loops. - The starting point is a dimcr 
representation of the wave function. As remarked by 
SutherlandcJ, the overlap < C\C > between two normal- 
ized valence-bond states C and C can be diagrammat- 
ically computed by constructing their transition graph. 
This graph is obtained by drawing on the same lattice 
the oriented bonds of C and C The result is a set of 
locally-oriented and non-intersecting loops which cover 
the lattice. Let n(C U C) be the number of loops and 
x(C U C) the number of bonds to be reversed for each 
loop to have alternating bond orientations. The result is 



< C\C' >= 2"( CuC "'~ A '/ 2 (— iy( CuC ') 



(6) 



As an illustration, consider the trivial case C — C . All 
loops are of length 2 and their number is n(CUC) = N/2. 
Each 2-site loop is already alternated and x{C U C) = 0. 
So we get < C\C >— 1. From the overcomplete family 
of valence-bond states one can formally extract a (non- 
orthogonal) basis. In this basis, the ground sate reads: 



|V>=]>>(C)|C 



> 



(7) 



2. No spin-Peierls symmetry breaking 

The N = 36 singlet sector of the J2 = —2 and J4 = 1 
Hamiltonian has a quasi degeneracy of the ground state 
(Fig. |l|). The lowest energy is a single state and the 
next one is 3-fold degenerate (impulsions k = Ai^,3). 
As mentioned before, the behavior of the singlet gap to 
the first Ai^.3 state (Fig. |l2|) suggests that these two 
levels are degenerate in the thermodynamic limit. The 
next pair of energies in S = may also be seen as asymp- 
totically degenerate k = 0, k = Aj states (Fig. p"3| ). 

The point is that no pattern gives a 4 or 8-fold degen- 
eracy compatible with the quantum numbers of these 4 
or 8 lowest eigenstates. In particular, any long-ranged 
ordered dimer covering leads at least to 12 degenerate 
states. For bigger patterns (4-site rhombus for instance) 
the number of broken symmetries is even bigger and so 
would be the degeneracy. Therefore, the 1+3 degeneracy 
of the two first eigenvalues (k = and k = ^1^,3) can- 
not be explained in such a conventional spin-Peierls like 
picture. 



The normalization of > gives 

< ^ >= £ [^(C) t ^(C')2" (CuC " ) (-l) x(CuC ' ) ] (8) 
c,c 

Small loops hypothesis. — This equation can be re- 
written in term of a sum over the transition graphs, or 
loop coverings Q of the lattice: 



< ^ >= 2 n{G \-i) x{ ^v{g) 



where we defined: 



c,c /cuc=e 



(9) 



(10) 



This is possible because any decomposition C U C of a 
loop covering Q has the same (— l^ x ( CuC ) = (_i) :c (&) a s 
well as the same n(C U C) — n{Q). We now turn to the 
important hypothesis: we assume that the only graphs 
Q which contribute to Eq. jl0[ ) are those which do not 
contain any loop surrounding the torus. This statement 
about the non-contribution of topologically non-trivial 
loops implies that there should not be any long distance 
coherence in the wave function. We stress that this is a 
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much stronger requirement that demanding short dimers 
(for instance, a Valence Bond Solid or a Valence Bond 
Crystal state does not fulfill this hypothesis) . In the MSE 
spin liquid all the correlations we measured (spin-spin, 
chiral-chiral and dimcr-dimer) seem to decay exponen- 
tially over a distance £ and we do nothing but assume 
that the loops contributing to Eq. (^) have also a finite 
characteristic length-scale. This strong condition implies 
the 4-fold degeneracy of the ground state. We propose 
the arguments summarized below. 

Cuts along the torus and J^-fold degeneracy. — The 
idea is now to chose a cut A surrounding the torus in 
one direction. Each valence bond-state C has a number 
A(C) of bonds crossing A and it is possible to change 
the sign of the configurations C for which A(C) is odd: 



W >=^(-i) A(c V(c)|c> 

c 



(11) 



If the sample has an odd number of sites in the direction 
perpendicular to A and if \tp > has an impulsion k, it 
is easy to check that the state \ip' > has an impulsion 
k' = k + A, where A is one middle of the boundary of 
the Brillouia zone (see inset of Fig. [l^). But we have 
also shownED that \ip' > has the same energy as \ip > 
provided \tp > only "contains" small loops (hypothesis 
above). This two- fold degeneracy must be extended to 
four when the sample has a 27r/3 rotation symmetry since 
irreducible representations for wave vector k = A are of 
dimeeeion 3. A similar construction was provided by 
WenE3 in the framework of a mean-field model for P— 
and T— symmetric spin liquid. In Wen's work, a state 
\ip' > is built from the original ground state > by 
a single gauge transformation on the bonds degrees of 
freedom. This transformation amounts to introduce a 
flux quantum through a non contractable loop A. 

Now we can interpret the N = 36 singlet sector in the 
following way: as the system becomes larger than £ the 
topologically non trivial loops grow scarce exponentially 
with y/~N/£. Consequently the lowest singlet eigenstate 
become 4-fold degenerate (k = and k = ^1,2,3). At 
N = 36 this topological degeneracy is still not very accu- 
rate but we can clearly distinguish two quasi quadruplets. 
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13. Lowest energy levels of the N = 36 sample, 
( J2 = —2 and J4 = 1) in the sectors of k = 0, A and B. The 
symbols represent the quantum number of the different eigen- 
states. For instance, the ground state (S — 0, E — —142.867) 
belongs to the trivial representation of the spatial group (full 
upward triangle, labeled k = 0, i?2x — l,R n = l,a — 1). 



4- A Chiral spin liquid ? 

Some time ago it has been suggested that an an- 
tiferromagnetic quantum spin liquid could be chiral 
and could ibjeak the time-reversal (T) and parity (P) 
symmetriesc3~E3. This scenario was definitely discarded 
in the case of thej-simple Heisenberg Hamiltonian on the 
triangular latticeEU. Th£L situation on the Kagome lat- 
tice is more ambiguousEj. In the MSE case, the pres- 
ence of a gap in the magnetic excitations and the de- 
generacy of the ground state in the thermodynamic limit 
would be consistent with a chiral theory. Indeed, WenU 
argue that an incompressible chiral fluid should have 
a 2/c 9 -fold vacuum degeneracy on a manifold of genius 
g (but this is not specific to chiral models: one also 
finds aA 5 topological degeneracy in P- and T-symmetric 
modelsc3). On the other hand, as underlined before, 
the MSE does not support LRO in the chiral parame- 
ter k — 2{S\ x S3) ■ S%. Therefore, the only possibility 
woulcLbe an hypothetical "Chiral Spin Liquid" ground 
statec3l23 where the expectation value of a cyclic permu- 
tation operator i , i ll i 2 ,i 3 , ...,»„ acquires a non-zero imagi- 
nary part on large loops i%, 13, . . . , i n . In such the- 
ories elementary excitations are unconfined spin-i ex- 
citations. The spectra of N = 27 and N — 28, are 
the largest pair of consecutive sizes we diagonalized and 
are believed to be close to the asymptotic limit. Their 
comparison does not plead in favor of this unconfined 
spinon hypothesis: the ground state energy per site of 
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N = 27 (E s=0 /N = -3.919) is slightly lower than the 
first triplet energy per site of N = 28 (E s=1 /N = -3.918 
and E s=0 /N — —3.963). This is rather a signature of ele- 
mentary spin-1 excitations and makes the chiral scenario 
unlikely. 

5. An AKLT / Haldane system ? 

We have shown, so far, that the MSE spin liquid is 
characterized by a low ground state degeneracy. For half- 
integer spins in one dimension, the LSM&3 theorem forces 
the ground state to be degenerate or gapless. Some at- 
tempts to generalize it in two dimensionsE-3 suggests that 
a ground state degeneracy could also be generic in a (half- 
integer) gapped state in two dimensions. From this point 
of view, our results in the MSE case are not unexpected. 
But such a degeneracy is usually associated to a sponta- 
neous breaking of a discrete symmetry, say translation, 
and a dimer or plaquette local order parameter. If our 
analysis of the N = 36 spectrum is correct, such a spin- 
Peierls phenomenon is absent. Without any evidence of 
a simple local order parameter, the topological point of 
view is natural and states that a ground state degeneracy 
can occur without any local order parameter breaking 
a discrete symmetry (a famous example is the Laugh- 
lin wave function, "phich discrete degeneracy depends on 
the space topologyE£l). However, the translation symme- 
try is almost certainly broken in the MSE case (at least 
four states with different impulsions are asymptotically 
degenerate) and one should be careful before excluding 
the possibility of any kind of LRO. In search for an al- 
ternative to the topological degeneracy, we looked for 
explicit wave functions (for half-integer spin) with a low 
degeneracy and the minimum number of broken spatial 
symmetries to compare with numerical results. 

For high spins matching the coordination number 
(2S = z) the AKLT wave function is non degenerate. 
When this criterion is not satisfied, nevertheless, the 
AKLT procedure can help to write wave functions with 
few broken symmetries. The idea is to chose a superlat- 
tice with n spins in the unit cell and of coordination num- 
ber z = n. On the superlattice, a unique AKLT state is 
written by identifying the n spins of a cell as a larger spin 
S = 5 . The state we obtained has a degeneracy which is 
the number of possibilities to define a superlattice L over 
the real (triangular) lattice. Consider the S = 2 AKLT 
wave function for the square lattice. It is mapped onto 
the spin— i triangular lattice by identifying the S = 2 
spins of the square lattice to 4 spin— ^ (on a rhombus) 
of the triangular one (Fig. |l4|.a). The choice of the po- 
sition of the square superlattice is done among twelve 
different possibilities: a configuration can be translated 
at four different places and rotated in three directions. 
One should remark that these states do have some long- 
ranged order. Let II a be the projection operator of the 
rhombus a in its S = 2 subspace. By construction, the 



8-point correlation function < n a IIb > is long-ranged. 

We computed the ground state energy inside this 12- 
dimensional AKLT subspace on a small system (N = 12 
at J4 = 1 and J2 = —2) and found a good variational 
energy. The exact ground state energy is E/N ~ —4.128 
and the AKLT energy is E/N = — 3.75. The energy scale 
is given by the energy of the ferromagnetic state {E/N — 
0), and the energy of the ground state E/N = —3 for clas- 
sical spins. This variational result is the best achieved so 
far and is even much better than the Schwinger-boson 
energy for this Hamiltonian (the Schwinger-boson so- 
lution is a two-sublattice (NLRO) collinear state with 
E/N = —2.96). The reason for this competitive en- 
ergy is that the AKLT construction naturally provides 
some local ferromagnetism inside a plaquette which low- 
ers J2-P1.2 and some antiferromagnetism at larger dis- 
tance (between neighboring rhombus) which insures a 
low J4 (Pi 2,3,4 + h.c). When the size of the unit cell is 
varied, one tunes the strength of the antiferromagnetic 
correlations at short distance. For instance, it is possi- 
ble to enhance it if we start from the wave function of 
the 5 = 3/2 honeycomb lattice (Fig. |l4|.b), where the 
unit cell has only 3 sites. In this case we find an en- 
ergy slightly higher (E/N = -3.42 on 12 sites). On the 
other hand, varying the ferromagnetic coupling J2 in the 
Hamiltonian could bring the ground state closer to an 
AKLT wave function. 

The twelve AKLT states are not coupled by any lo- 
cal Hamiltonian in the thermodynamic limit, and should 
therefore be degenerate in an infinite system. In a theory 
where all excitations are gapped, the low energy physics 
(except for possible edge states) is determined by the 
degenerate ground state. Therefore, this twelve-fold de- 
generacy characterizes the fixed point and we expect the 
same degeneracy in the MSE problem if it belongs to this 
universality class. Unfortunately the quantum numbers 
only partially coincide with the lowest eigenstates found 
in the N = 36 spectrum. So, despite of a good varia- 
tional energy on N — 12, these particular trial states do 
not provide a quantitative explanation of the 4 or 8-fold 
degeneracy observed on 36 sites. However, the family of 
these extended AKLT states certainly captures a part of 
the local constraints imposed by frustrating MSE cou- 
plings. Furthermore we cannot completely exclude that 
the N = 36 might not be a faithful picture of the infinite 
system spectrum and that the ground state could have 
a complex Valence Bond Solid order of the type we have 
just described. 
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FIG. 14. Left: black rhombus represent the S = 2 spins of 
the AKLT wave function. They form a square lattice (dashed 
lines). Right: honeycomb superlattice of 3-site plaquettes on 
the triangular lattice. There are 18 different ways to place this 
honeycomb superlattice on the triangular one (3x2 transla- 
tions and 3 rotations). 



V. SUMMARY AND CONCLUSIONS 

We analyzed the spin liquid phase of MSE model on the 
triangular lattice. Due to important frustration, LRO is 
destroyed by quantum fluctuations at zero temperature. 
The system has short-ranged correlations: < Si ■ Sj >, 
dimer-dimer and chiral-chiral correlations probably de- 
cay exponentially with distance. The spectrum has a 
spin gap and the comparison between odd and even sam- 
ple pleads in favor of spin-1 excitations. As for the spatial 
symmetries and degeneracy of the ground state (presum- 
ably 4 or 8), they provide evidence of no spin-Peierls 
nor simple plaquette order. This is a generic feature in 
systems where the spin matches the coordination num- 
ber (25* = z) and where an AKLT wave function can 
be constructed, but it is unconventional for spin-i. We 
proposed two interpretations. The first possibility is an 
AKLT or Haldane -like phase. We constructed such a 
wave function, which is not a tensor product of plaquette 
states. Present numerical data on the ground state sym- 
metries do not completely agree with this particular trial 
state but, because of its very low variational energy, such 
a scenario cannot be excluded. Secondly, this degeneracy 
might be a consequence of the non trivial topology im- 
plied by periodic boundary conditions. We emphasized 
that it must be the case if the system has absolutely no 
long-ranged coherence. 

This spin liquid phase might be the very first magnet 
with a "disordered" and gapped quantum ground state 
with no simple local order parameter in two dimensions. 
From this point of view it can be seen as a prototype of 
the RVB state as was proposed by AndersontZl. 
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